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Abstract 

We study the properties of a new class of gamma-ray sources consisting of millisecond 
pulsars totally or partially surrounded by evaporating material from irradiated companion 
stars. Hidden millisecond pulsars offer a unique possibility to study gamma-ray, optical 
and radio emission from vaporizing binaries. The relevance of this class of binaries for 
GRO observations and interpretation of COS-B data is emphasized. 

1. Introduction 

The missing link between low mass binaries (LMXBs) containing rapidly rotating neu- 
tron stars and binary millisecond pulsars has been provided by the discovery of pulsars 
in the process of vaporizing their very low-mass companion star. Two such systems are 
currently known: PSR. 1957+20 with spin period P — 1.6 msec (Fruchter, Stinebring and 
Taylor, 1988) and PSR 1744-24 A with P = 11.6 msec in the globular cluster Terzan 5 
(Lyne et ah, 1990; Nice et ai, 1990). PSR 1957+20 is the prototype of SVPs with stable 
eclipse properties (type I) whereas PSR 1744-24A shows a more erratic and sometimes 
intermittent behavior of radio-eclipses characteristics of a second ciass of SVPs (type II). 
In addition to these two SVP-classes, a third class of SVPs is expected on theoretical 
grounds (type III). The SVPs of type III are completely enshrouded by the evaporat- 
ed material from the irradiated companion star, and they are therefore ‘hidden’ pulsars 
because their radio-emission may be completely blocked. However, they might be re- 
vealed by X-ray /gamma- ray emission produced by the interaction of the pulsar wind with 
evaporating material at the termination shock (Tavani, 1991b, c). 

The existence of ‘star- vaporizing pulsars’ (SVPs) is in agreement with an evolutionary 
scenario of low-mass binaries containing compact stars whose energetic radiation (X-rays 
and 7 -rays) illuminates their companion stars and drives a strong evaporative mass out- 
flow during the LMXB and post-accretion evolutionary phases (Ruderman, Shaham and 
Tavani, 1989, hereafter RST; Tavani, 1989, 1991a, b,c). Further ‘vaporization’ of the com- 
panion star in the post-accretion phase is caused by radiation produced by the interaction 
of the pulsar wind with the outflowing material from the companion. It is precisely this 
evolutionary phase that leads to the formation of a new class of 7 -ray sources consisting 
of hidden millisecond pulsars. We follow here the discussion of hidden millisecond pulsars 
presented in Tavani (1991c), with emphasis on their gamma-ray emission (Tavani, 1991f). 
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2. Three Classes of Star- Vaporizing Pulsars 

Pulsar emission in the form of a Poynting flux and a relativistic particle wind is 
possible if r\ c < 77, with r/ c the light-cylinder radius r/ c — c/fi = (4.77 • 10 6 cm) * P„ 3 , 
where P_ 3 = P/(10” 3 sec), and 77 the distance from the pulsar at which the material 
outflowing from the companion star is stopped by radiation or magnetic pressure. For the 
cases considered in this paper, 77 = a — with a the orbital distance, and R s the distance 
from the companion star at which the ram pressure of the outflowing gas is balanced by the 
radiation pressure from the pulsar. For SVPs, r $ is in the range r/ c < r s < a, and satisfies 
the relation I//(47rr 2 c) = amv g l{2ixr l 2 ) with L the pulsar luminosity, r the distance from 
the pulsar, c the speed of light, m = |m| rac i the mass loss from the irradiated companion 
star, r f — a — r 3 = R s the distance from the companion, a a dimensionless quantity giving 
the hydrodynamical density enhancement near R s relative to the density obtained from 
a \/r n ideal dependence for constant wind velocity 1 , and v g the gas velocity. Since the 
wind velocity of radiation-driven mass outflows in SVPs and LMXBs is of the same order 
of magnitude of the escape velocity from the star v e (Tavani and London, 1991) and since 
the Keplerian orbital velocities vk of low-mass binaries with orbital periods P 0 r & ~ 1 day 
are comparable with u e , we can write the pressure balance equation as follows 

fCF[x] — m (1) 

where we defined x = r 5 /a, fC = L yja / ( 2ca y/ GM \ ) ~ (1.53T0 17 g s~ 1 ) y /a n P^ 2 o ^45/(<*^ 3 )> 
with jP_2o given by P = 10 _2 ° P-20 s 5” 1 , I4 5 the neutron star moment of inertia in units 
of 10 45 g cm 2 , a = 10 11 an cm, and Mi the pulsar mass. The function F[x ] is given by 

F[x] = (1-.t) 3 / 2 x- 2 (2) 

Eqs. (1) and (2) give the position of the shock radius as a function of rn , orbital (an, M \ ), 
pulsar (P- 20 , P-s) ^ n d hydrodynamics (a) parameters. 

It is also useful to write the radiation-driven mass loss rate (in cgs units) as (RST) 
m 10~ 17 f L A = (3.95 * 10 16 g s~ l ) f P-20I45A-2/ P-3^ where the quantity / depends 
on the net heating and hydrodynamics of the mass outflow near the companion star 
(Tavani 1989), and with A = [R^/(2a)] 2 (and A = 10 _2 A_2) the solid angle subtended by 
the companion star to the primary where R+ is the companion radius. For typical fluxes 
irradiating SVP companions we have / ~ 0.1T, where T is a dimensionless spectral 
‘quality’ parameter which gives the efficiency of energy deposition. For the two energy 
bands which are effective in driving a strong wind, we obtain T ~ 10 — 100 for soft X-rays 
and T ~ 0.2 for hard X-rays and 7-rays (Tavani, 1989). 

We can then rewrite the pressure balance condition Eq. (1) in a form which, in first 
approximation, is independent of the pulsar luminosity and depends upon the spectral 
‘quality’ of the heating / as well as upon a and orbital parameters 

F[x] = 0.25S/aA_ 2V /ji77 a" 1/2 (3) 

1 Results of a hydrodynamical calculation of the outflow in the case of PSR 1957-j-20 and PSR 1744- 
20A gives 1 ~ a ~ 10 (Tavani and Brookshaw, 1991). 
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Given the variety of pulsar luminosities and orbital parameters of low-mass binaries we 
can expect several types of SVPs characterized by different R,'s and shapes of the mass 
outflow influenced by the combined action of gravitational, radiation and Coriolis forces. 

SVPs of type I are obtained when m is close to mi = /CF’fl — R,/a] and the Coriolis 
forces acting on the wind particles near R s are not large enough to produce an appreciable 
tangential component of the gas pressure. The radio eclipses expected for type I SVPs are 
sharply defined and possibly symmetric (see Fig. 1). In the case of PSR 1957+20 (the 
prototype of type I SVPs) we obtain m x ~ 10 14 g s _1 for Rja = 0.0353. A lower limit 
for / in the case of PSR 1957+20 is obtained by equating m with rhi, i.e., / ~ 0.213/a. 
The upper limit on / is obtained by the relation / ~ (0.213/a)(F 1 [l — R s / a]/ F^l — R,/a]) ~ 
10/a, where we used the shock distance from the companion appropriate for PSR 1957+20 
to be R 3 /a ~ 1/5 as suggested by hydrodynamical calculations of the mass outflow fitting 
the eclipse characteristics of PSR 1957+20 (Tavani and Brookshaw, 1991). 



FIG. 1 - Schematic behavior of SVPs of different types (based on the hydrodynamical calculation 
of Tavani and Brookshaw (1991)). The asterisk denotes the position of the pulsar and the 
companion star is assumed to be rotating counterclockwise. 
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SVP’s of type II and III are possible when m is larger than rhi and the Coriolis 
forces together with outflow velocities of order of Keplerian velocities are sufficiently large 
to produce a large tangential component of the gas pressure near R s . This additional 
component of the gas pressure forces the gas to begin to surround the pulsar in the 
leading edge of the eclipsing region. The eclipses of type II SVPs are expected to be 
occasionally non-symmetric with a time variable leading edge of the eclipse (see Fig. 1). 

A remarkable example of SVP of type II is offered by PSR 1744-20A which shows 
symmetric (Nice et a/., 1990), non-symmetric (Lyne et al,, 1990) and total eclipses. We 
obtain in this case mi ~ (3.71 • 10 12 g s~ l )P_ 2 o$ 2 / a > where 4/ — RjRi, ~ 1 with Ri the 
Roche lobe radius, and R*/a 0.18$ (Lyne et al, 1990). As in the case of PSR 1957+20 
we obtain for PSR 1744-20A the possible range for /, 0.16/a ~ ~ 10/a, where the 

upper limit has been obtained by taking Rs/a ~ 1/2, as suggested by the hydrodynamical 
calculation of the outflow of PSR 1744-20A (Tavani and Brookshaw, 1991) . 

3. Hidden Millisecond Pulsars 

SVPs of type III are obtained if the material near R s has sufficient tangential kinetic 
energy to overcome the pulsar radiation pressure and is able to surround the pulsar 
completely (see Fig. 1). Eq. (3) shows that this case may occur for a favorable combination 
of the quantities /, a, A, and a. PSR 1744-20 A offers a good example of a SVP of type 
III during its ‘radio-quiescent’ mode observed several times in 1990 (Lyne et al , 1990; 
Nice et al, 1990). PSR 1744-20A is an example of a SVP which is able to switch 
between different eclipse modes, but SVPs which spend most of their time as type III 
SVPs may exist as well. In this case the central pulsar is completely surrounded by the 
wind material in the orbital plane that partly accumulates near R s and partly escapes 
from the binary. This class of ‘hidden’ millisecond pulsars may have their radio emission 
completely quenched unless the wind material is confined in a thin disk. The shape of the 
enshrouding material in three dimensions depends on the interplay among gravitational, 
radiation and pressure forces, and a ‘bubble’ or a donut-shaped ‘torus’ can describe the 
gas configuration surrounding the pulsar. The hydrodynamics together with the stability 
and time dependent effects influencing the enshrouding material is the subject of intense 
research (Tavani, 1991c,f) and an account of the results is beyond the scope of this paper. 
Furthermore, even though the currently known SVPs and the related numerical modelling 
suggest the existence of a gas configuration enshrouding the pulsar which is optically 
thin to hard X-rays and gamma-rays (Tavani and Brookshaw, 1991), a more complex 
configuration may be possible. In the following we will address the main qualitative 
features concerning the emission of gamma-rays expected from hidden millisecond pulsars 
assuming that the evaporating material is optically thin to hard X-rays and 7-rays. 

4. Gamma rays from shock acceleration 

Rapidly rotating pulsars are believed to produce a radially expanding relativistic wind 
of ions and e^pairs in addition to the comoving Poynting flux of electromagnetic energy 
in a MHD approximation (e.g., Kennel and Coroniti, 1984; Hoshino et al , 1991). The 



MHD pulsar wind is characterized by the Lorentz factor of particles and by the ratio a of 
electromagnetic to kinetic energy of wind particles. We can express the Lorentz gamma 
factor 71 in terms of a parameter 77 which gives the efficiency of acceleration and pair 
production near the pulsar as 

71 = 7](Ze$ c )/(miC 2 ) ~ 2.3 • l0 6 r,(Z/A)(P. 2O / Pl 3 ) 1/2 . (4) 

where $ c ~ [2ty R p j cP) 2 {2 /.tj R^)R P is the voltage across open field lines of a pulsar of 
magnetic moment //, and Z/A is the average ratio of charge over mass number of the 
wind particles. The quantity a is 

a = B 2 /4nc 2 [Nirrii'fi + (iV e + + N e -)m e 7 e ±] (5) 

where TV,-, m t - , 7,- and N e ± , m e , 7 e ± are the number densities, masses and relativistic Lorentz 
factors for the ions and e^pairs, respectively. In the MHD approximation, 7 e ± = 7 = 71 . 
Detailed models for the Crab nebula require a ~ 10 -2 — 10 -3 (Kennel and Coroniti, 
1984), and 7/ ~ 0.3 (Hoshino et a/., 1991) which implies a value of 71 between 10 5,3 (fully 
stripped iron) and 10 5 * 6 (protons). 

The pulsar wind of relativistic particles and Poynting flux interacts with the inner 
boundary of the SVP ‘bubble 1 or ‘torus 1 of Fig. 1 and a relativistic shock ensues. For 
simplicity, we will assume that the region around the shock region is sufficiently dense 
that the e.m. waves and plasma in the pulsar wind do not penetrate the mass outflow, and 
that the magnetic field is transverse to the shock normal at r s . Under these conditions, 
diffusive Fermi acceleration (e.g., Ellison et a/., 1990) is not applicable, and the shock 
is transverse. The characteristics of particle acceleration depend on the details of the 
relativistic shock and several models are currently being investigated (Tavani, 1 99 If). 

One sufficiently developed model of shock acceleration based on resonant absorp- 
tion of magnetosonic waves by positrons near the ion reflecting region has been re- 
cently applied to the study of gamma-ray emission from PSR 1957+20 (Arons and 
Tavani, 1991). The radiation spectrum emitted downstream depends on the physical 
processes regulating the efficiency of synchrotron, inverse Compton scattering and non- 
thermal emission of the relativistic shock. The typical synchrotron energy is of order 
E 3yn ~ (10 keV)'yl ± 6 Bp y Ja/(l + a), where 7 e ± = 10 6 7 e ±,6 and Bp is equal to the mag- 
netic field of an ideal Poynting flux assumed to carry all the spindown energy of the pulsar. 
Typically, Bp = fi/(rf c r) ~ (440 G)B 9 /(P 2 3h r, tll ), with B 9 the surface magnetic field in 
units of 10 9 G, and r Srll = r 5 /(10 n cm). 

The radiation spectrum expected from the relativistic shock with ion-induced acceler- 
ation of positrons has been shown to have a photon index s ~ 2 (Hoshino et a/., 1991). 
The effective synchrotron emissivity of the shock is a power law j^/e oc e“ 2 ( in units of 
photons sec” 1 cm -3 ), extending from to e[ b \ By equating the radiation loss timescale 
with the ion-induced acceleration timescale we obtain a range for the upper energy cutoff 
10 MeV ~ e[ b ) ~ 10 GeV for different choices of 7! and a. As an example, an application 
to PSR 1957+20 of the ion-induced acceleration model with pulsar wind parameters for 
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the Crab nebula gives t [ b) ~ 3 MeV (Arons and Tavani, 1991). If ions contribute to a 
substantial part of the pulsar wind energy, the radiation spectrum is expected to be a 
power law of index s from a few keV to tens or hundreds of GeV. The corresponding 
heating on the surface of the irradiated companion yields 10~ 3 ~ / ~ 10 -2 . Alternately, 
if ions do not appreciably contribute to the pulsar wind kinetic energy and therefore are 
not effective in pair acceleration near r s , most of the pulsar wind energy is expected to 
be radiated at lower energies peaked around E syn (Hoshino et al., 1991) unless an ad- 
ditional shock acceleration mechanism (such as shock-drift acceleration) operates. The 
corresponding spectrum is expected to be in the X-ray energy range and 1 ~ / ~ 0.1. For 
energies below e 1? the spectrum is cut off, since the emission comes from the low energy 
tail of the nonthermal part of the particle spectrum, and from the Maxwellian pairs at 
energies less than 71 nr±c 2 , with tj ( oc e 4 ^ 3 (e.g., Jones and Hardee, 1979). 


5. Observational Consequences 

Hidden millisecond pulsars produce a gamma-ray yield larger than in the case of 
isolated pulsars, and they are therefore interesting for interpreting galactic gamma-ray 
sources. The 7-ray flux from hidden millisecond pulsars at the distance D can be estimated 
for a power-law emission as follows (t is the photon energy) 


dN{E > e) 
d t 


(3-10 4 ph. cm 2 s 4 ) 


■‘p, 36 


i 1 MeV 


Dl P c 0-2 


( 6 ) 


where £ is the efficiency of conversion of pulsar wind energy into radiation (Hoshino et 
al ., 1991), L Pi 36 = L p /{ 10 3r> erg s -1 ), and D kpc = D/( 1 kpc). 

The detectability of gamma-rays from a hidden millisecond pulsar depends crucially 
on the ratio L Pi 36 / D\ . In the case of GRO detectors, we notice that hidden pulsars 
would be detectable by OSSE, COMPTEL and EGRET for a normal pointing time of 
two weeks if L p ^ & ! D \ pc were larger or equal to ~ 1, 0.2, and 0.01, respectively. 


5.1 Hidden Pulsars with Main Sequence Companions 

Low-mass X-ray binaries (LMXBs) are believed to be the most probable progenitor of 
binary millisecond pulsars (e.g., Tavani, 1991a). However, no galactic binary millisecond 
pulsar has yet been discovered with a main sequence companion (e.g., Tavani, 1991e) 
even though the majority of LMXBs have orbital periods proper to binaries with main 
sequence companions (e.g., Ritter, 1990). This might be due to the effect of Doppler 
spread of the radio signal in orbits with orbital periods P or b ~ 6 hrs. However, millisecond 
pulsars with main sequence companions satisfy Eq. (3) with A_ 2 ~ 1, and they may 
be hidden SVPs. Hidden SVPs with main sequence companions may be revealed by 
the X-ray /7-ray emission from the radio-obscuring ‘bubble’ and, possibly, by an optical 
modulation of the irradiated companion star if the inclination angle is favorable. The 
luminosity of the companion star reprocessing the secondary radiation from the pulsar 
wind is L ~ (4 • 10 33 erg s -1 )P_ 2 ol45 A_ 2 /P 3 3 with a photospheric temperature T ~ 
(1.5 • 10 4 ° K)P- 2 qP- 3 ^ 4 a n /2 - Main sequence companion stars of millisecond pulsars 
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would appear brighter and ‘bluer’ compared to similar undisturbed stars. 

5.2 Interpretation of COS-B sources 

Even though the pulsed and unpulsed radio emission of hidden SVPs is probably 
completely blocked, they can be visible in the X-ray/ 7 -ray band and, indirectly, in the 
optical band because of irradiation of the companion star (see Sect. 5.1). The radiation 
spectrum originating from the ‘bubble’ surrounding a SVP can be a power-law extending 
to energies ~ 10 GeV, and hidden SVPs can be important for the interpretation of X-ray 
and 7 -ray sources already discovered by COS-B and/or to be discovered by the Gamma 
Ray Observatory (GRO). The association of an X-ray and/or 7 -ray source with a ‘blue’ 
companion star possibly showing optical modulation provides a well defined signature 
for the existence of hidden SVPs (Tavani, 1991c). The detectability of this association 
depends on details such as the inclination angle, and the amount of shadowing of the 
optical emission. It is however worthwhile to search for optical modulations in the field 
of known 7 -ray sources (e.g., unidentified COS-B sources , Swanenburg et a/., 1981). 

Furthermore, again depending on the geometry, even the gamma-rays can be modu- 
lated with an orbital period by scattering and absorbing material related to either the 
companion star itself or to density enhancements in the bubble/ torus surrounding the 
pulsar. A well defined signature for a successful association of gamma-ray sources with bi- 
nary systems is therefore the existence of modulations in the gamma-ray flux with periods 
of few hours. A search of orbital period modulations in the gamma-rays from COS-B-like 
sources has not been systematically attempted before, and a periodicity search of GRO 
and COS-B data is strongly urged. Folding gamma-ray data with given periodic patterns 
may significantly increase the signal to noise ratio. 

5.3 Gamma rays from globular clusters 

A relatively large number of single and binary millisecond pulsars (~ 25) has been 
discovered in globular clusters (e.g., Tavani, 1991e), possibly as a consequence of the 
high probability of stellar encounters in dense cluster cores. A remarkable example of a 
globular cluster which harbors many millisecond pulsars (~ 10) is 47 Tucanae (Manchester 
et a/., 1991). It is therefore natural to expect that a relatively large number of pulsars 
in globular clusters are hidden, and that a detectable gamma-ray or optical signal can 
be used to deduce their existence. We note here that a class of ‘blue stragglers’ as those 
discovered in NGC 6397 (Auriere, Ortolani and Lauzeral, 1990)and in 47 Tucanae (Paresce 
et a/., 1991) might contain millisecond pulsars whose irradiation alter the properties of 
their companions (Tavani, 1991b, d). The existence of pulsar-driven blue stragglers is 
particularly interesting in globular clusters where the majority of stars are main sequence 
(Tavani, 1991d). 

We note that the gamma-ray flux from individual and not enshrouded pulsars in 
globular clusters is too low to be detected by current gamma-ray detectors. However, 
the shock acceleration mechanism of hidden powerful pulsars enhances the gamma-ray 
yield for each single source, and single (in the case of EGRET) or collective emission from 
hidden pulsars might be detected by GRO. A periodicity search in a range of a few hours 
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might enhance the discovery probability. We estimated the number of sources necessary 
to make globular clusters detectable by the GRO instruments. We assumed the photon 
flux from Eq. (6) and an average Lp , 36 = 1. For 47 Tucanae we obtain the minimum 
numbers (21, 4 , 1) for OSSE, COMPEL and EGRET, respectively. For Terzan 5 we 
obtain (50, 10, 3), and for the core collapsed cluster NGC 6397 we estimate the minimum 
numbers in the range (5, 1, 1). 

Work was performed at IGPP-LLNL under the auspices of the U.S. DOE under con- 
tract W-7405-ENG-48 with NASA support, grants GRO-90-29 and GRO-PFP-91-23. 
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